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ABSTRACT 

We use the spherical collapse (SC) approximation to derive expressions for the 
smoothed redshift-space probability distribution function (PDF), as well as the p- 
order hierarchical amplitudes S p , in both real and redshift space. We compare our 
results with numerical simulations, focusing on the Q = 1 standard CDM model, 
where redshift distortions are strongest. We find good agreement between the SC 
predictions and the numerical PDF in real space even for ox > 1, where <jl is the 
linearly-evolved rms fluctuation on the smoothing scale. In redshift space, reasonable 
agreement is possible only for ol < 0.4. Numerical simulations also yield a simple 
empirical relation between the real-space PDF and redshift-space PDF: we find that 
for er < 1, the redshift space PDF, P [$( z )]> is, to a good approximation, a simple 
rescaling of the real space PDF, P [5], i.e., P [S/a] d[5/a] = P w 2 )Mj) d [<5(z)/c(z)] , 
where a and ov z ) are the real-space and redshift-space rms fluctuations, respectively. 
This result applies well beyond the validity of linear perturbation theory, and it is a 
good fit for both the standard CDM model and the ACDM model. It breaks down for 
SCDM at ex ~ 1 , but provides a good fit to the ACDM models for a as large as 0.8. 
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1 INTRODUCTION 

One of the central statistical quantities of interest in the 
study of large-scale structure is P[p], the one-point probabil- 
ity distribution function (PDF) of the density field p. In the 
linear regime, the PDF simply retains its Gaussian shape, 
with the variance scaling up as the square of the growth 
factor. In the nonlinear regime, however, P[p] diverges away 
from its (presumably Gaussian) initial shape. 

A number of recent studies have been undertaken to 
predict analytically the evolution of the PDF. Kofman et 
al. (1994) calculated the PDF produced by the Zeldovich 
approximation. Given a formalism to calculate the cumu- 
lants of the evolved distribution (Peebles 1980; Fry 1984; 
Bernardeau 1992a, hereafter B92), the Edgeworth expansion 
can be employed to derive an approximate expression for the 
PDF (Juszkiewicz et al. 1995, Bernardeau & Kofman 1995). 
A better alternative is to derive a similar expansion around 
the Gamma PDF (Gaztanaga, Fosalba & Elizalde 2000). Un- 
like the Edgeworth expansion, the Gamma expansion always 
yields positive densities. The generating function for these 
cumulants can be used to derive the Fourier transform of 



the PDF, which can be inverted to give the PDF in the form 
of an integral (Bernardeau 1994). Alternatively, this gener- 
ating function can be transformed into a local Lagrangian 
mapping which, when applied to the initial Gaussian den- 
sity field, gives a PDF with the desired tree-level hierarchical 
amplitudes (Protogeros & Scherrer 1997, hereafter PS; Pro- 
togeros, Melott, & Scherrer 1997). Fosalba and Gaztanaga 
(1998a, hereafter FG) showed that when shear is neglected, 
the equations of motion automatically produce a local La- 
grangian mapping which exactly reproduces the tree-level 
hierarchical amplitudes (also Fosalba and Gaztanaga 1998b, 
Gaztanaga & Fosalba 1998). It is this approach, neglecting 
shear in the cosmological equations of motion, which we will 
utilize here. Fosalba & Gaztanaga dubbed their approxima- 
tion the "spherical collapse" (SC) approximation; although 
we adopt their terminology, it would perhaps be more accu- 
rate to refer to this approximation as "shear- free collapse" . 
The SC model has also been used to model the PDF of QSO 
Ly-a absorbers (Gaztanaga & Croft 1999). 

Of course, the three-dimensional distribution of galaxies 
is observed in redshift space rather than real space, with 
the radial coordinate distorted by the peculiar velocities of 
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the galaxies. Hence, the most physically relevant quantity 
is not the PDF in real space, but the redshift-space PDF. 
While the evolution of the hierarchical amplitudes in redshift 
space has been examined in some detail (Bouchet et al. 1995; 
Hivon et al. 1995; Scoccimarro et al. 1999), comparatively 
less work has been done on the distortion of the full PDF in 
redshift space. Kaiser (1987) showed that, to lowest order, 
the rms fluctuation in redshift space, <J( Z ), is related to the 
real-space rms fluctuation a, via 



2 



I+3/0 



5 /n 



(1) 



where /n ~ Q 6 . (Eq. 1 assumes that the bias parameter is 
unity. We make this assumption throughout, although bias 
can be incorporated into our formalism if the form of the 
bias is prescribed). The effect of redshift distortions on the 
skewness was calculated by Bouchet et al. (1995) and Hivon 
et al. (1995). They found that redshift distortions have little 
effect on the skewness. For fi = 1, for instance, the redshift- 
space skewness is well-approximated by (Hivon et al. 1995) 
5*3( z ) = 35.2/7 — 1.15(n + 3); this compares with a real- 
space result of 5*3 = 34/7 — (n + 3). (This theoretical result 
for the redshift-space skewness does not agree with numer- 
ical simulations except at very large length scales; see, e.g., 
Fig. 3 below). A systematic method to calculate the higher- 
order perturbation theory kernels in redshift space has been 
given by Scoccimarro, Couchman, & Frieman (1999). The 
most detailed treatments thus far have been given by Hui, 
Kofman, and Shandarin (2000), who calculated the redshift- 
space PDF in the context of the Zeldovich approximation, 
and by Watts and Taylor (2000) who calculated the redshift- 
space PDF in second-order perturbation theory. 

In the next section, we derive the local Lagrangian map- 
ping corresponding to the SC model, and apply it to Gaus- 
sian initial conditions to derive the real-space density distri- 
bution function. In Section 3, we derive an approximation 
to the redshift distortion of the density field based on the 
SC approximation. We calculate the corresponding redshift- 
space hierarchical moments and the redshift-space PDF in 
the SC approximation and compare the PDF to numerical 
simulations in redshift space. We also derive an interest- 
ing empirical relation between the real-space PDF and the 
redshift-space PDF. Our conclusions are summarized in Sec- 
tion 4. 



2 THE PDF IN THE SC MODEL 
2.1 Theory of the SC PDF 

One class of approximations to the evolution of large scale 
structure are local Lagrangian mappings, in which the den- 
sity at a Lagrangian point q at a time t is taken to be a 
function only of t and the initial value of the density at q, 



!Z(q,*)=JV(t)/[<S £ (q)], 



(2) 



where r\ = p/p, with p the mean density, and 8l is the 
linearly-evolved initial density fluctuation (5l — rjL — 1): 



5 L (t) = D(t)S(to) 



(3) 



In equation (g), N(t) is a normalizing factor chosen so that 
/ P[r,]d V = 1: 



N(t) = 



f(*L 



(4) 



A stochastic field of density fluctuations <5(x) can be 
described in terms of its cumulants n p , which are given by 
k 2 = a 2 = (S 2 ), k 3 = (S 3 ), k 4 = <<5 4 ) - (S 2 ) 2 - An initially 
Gaussian density field evolving gravitationally has the prop- 
erty that K P /(T 2( ' P_1 ' goes to a constant in the limit where a 
goes to zero, where a is the rms density fluctuation (Peebles 
1980, Fry 1984, Bernardeau 1992). Hence, it is conventional 
to define the hierarchical amplitudes S p given by: 



S p {a) = 



(a 2 )p~ 



(5) 



The values of S p in the limit where a 



are called the 
tree-level hierarchical amplitudes. 

This formalism can be exploited to derive an approx- 
imation to the PDF in the quasi-linear regime. For small 
fluctuations, a generic local Lagrangian mapping can be ex- 
pressed in a Taylor series: 



*?(q) = /[&(«!)] 



(6) 



It can be shown (PS, FG) that for the local Lagrangian 
mapping given in equation expanded as a Taylor series 
as in equation (m), the hierarchical amplitudes are given by 



5 3 (0) = 3V2, 
S 4 (0) = 4u s + Ylv\ 



(7) 



and so on. These amplitudes can be made identical to those 
of exact tree-level perturbation theory if we take the map- 
ping in Eq. (|) to be given by (B92, FG): 



N 



9 (6>-sin<9) 2 
2 (1 - cos6») E 



Sl = |[|((9 -sine)] 2 / 3 , 
for 5l > 0, r\ > 1, and 

r 9 (sinhfl- Of 



(9) 



n 



N 



2 (cosh 
3 r 3 



5 l 4 



(slab. 0-0)] 



2/3 



(10) 



for 8l < 0, T) < 1. [In this paper, we assume Gaussian initial 
conditions, but this result can be extended to non-Gaussian 
initial conditions as well (Gaztanaga & Fosalba 1998)]. 

There is a physical basis for this result: in the limit 
where shear is neglected, the equations of perturbation 
growth become (FG) 

S + -5-^j 1 ^=4nGa 2 5(l + 5), (11) 
a 3(l + o) 

where the dot denotes the convective derivative with respect 
to conformal time. (Hence, equation (hlh is an equation for 
the Lagrangian density). But equation (pi) is identical to the 
equation for S in spherical collapse, so the solution is the lo- 
cal Lagrangian SC mapping given by Eqs. (^)- ( |lo| ) . Further- 
more, the Fourier-space integrals which give the tree-level 
hierarchical amplitudes are spherically symmetric, guaran- 
teeing that the SC approximation gives the correct hierar- 
chical amplitudes at tree level (FG) . Although this approxi- 
mation was dubbed the "spherical-collapse approximation" 
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(FG), it does not presume or imply that the actual physical 
collapse is spherically symmetric; rather, it proceeds from 
the approximation that the clustering is shear-free. This ap- 
proximation should not be confused with the well-known 
analytic results for actual spherical collapse (Peebles 1980). 
For the SC approximation we have: 



™ ~ 1.62, 
21 



682 

189 
446440 
43659 " 
8546480 
243243 



3.61, 



10.22, 
- 35.13, 



(12) 



and note that V\ = 1 to produce the linear theory result. 
As noted, these numbers give the correct contribution to 
the pth-order cumulants, k p , for Gaussian initial conditions, 
e.g., S 3 = 34/7 (Peebles 1980, Section 42). 

We now derive the PDF in the SC approximation. The 
SC local Lagrangian mapping which relates 5i(q) to 7?(q) 
is given by Eqs. (9) and (|lo|). In our discussions below, the 
mapping r\ — NJ(Sl) is always understood to refer to equa- 
tions (^) and (|Io|), although one cannot write down a non- 
parametric form for /(<5l). These equations also define an 
inverse function: Sl = fZ^iv/N); it is 5h{rj/N) as defined 



by equations M an 



which actually enters into our 



equations for the PDF. The evolved PDF (given, as usual, 
in Eulerian space) is (PS) 



(13) 



where Po is the initial (unevolved) density distribution func- 
tion. The normalizing factor TV is most easily derived from: 



N 



Pa{8 L {x)) j^-dx, 

x ax 



(14) 



where we have used the change of variables x = T]/N. For 
Gaussian initial conditions, Pq is a Gaussian, and we get 



P[rj] = 



2-KCJL 



exp 



S L (y/Nf 
2o\ 



1 d5 L {r)/N) 
rj drj 



where ox is the linearly evolved rms fluctuation: 
a L (t) =D(t)a , 
and N is given by 



TV = 



2-KIJL 



exp 



Sl(x) 2 



2al 



1 dS L (x) 
x dx 



dx. 



(15) 



(16) 



(17) 



This expression gives the PDF for an unsmoothed 
density field, which is not physically observable. If we 
smooth the final density field with a spherical top-hat win- 
dow function, then we obtain a new "smoothed" local La- 
grangian mapping density /s(5l), which is given in terms 
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of the unsmoothed mapping / through the implicit relation 
(Bernardeau 1994; PS; FG): 



r (z) = r + MzX. 



(24) 



fs(6x.) = f[ 5 L 



<T( fs(5 L )Ro) 



(18) 



where Ro is the radius of the spherical top-hat, and a(Ro) is 
the rms fluctuation of the linear density field (in practice, de- 
rived by integrating over the linear power spectrum) . Taking 
rjs = N/s(Sl), where N is the new (smoothed) normaliza- 
tion factor, this expression can be recast in the form 



5l(vs/N) =f- 1 (rjs/N) 



a(Ra) 



a(R QVs /N) ' 



(19) 



If the linear power spectrum is a power law, P(k) oc k n , 
equations (|lq) an d (|l£j) simplify to (Bernardeau 1994) 



-(n+3)/6] 



fs(S L ) = f [S L f s (S L ) 
which gives 

S L ( VS /N) = r\vs/N)[ Vs /N]^ +3)/6 



(20) 



(21) 



The expression for Sl from either equation ( 13) o r equation 
( pl| ) can be substituted into equations jl5[)-(|l7|) to obtain 
Ps(v) (the PDF of the smoothed density field). 



2.2 Test with Nbody simulations 

To test our approximations for the PDF, we compare to a set 
of 10 standard cold dark matter (SCDM) simulations with 
64 3 particles over an L — 180 Mpc/ft box. (We also used 
simulations with an L ~ 400 Mpc//i box to test for volume 
effects. See Baugh, Gaztanaga & Efstathiou 1995, for more 
details on these simulations). By smoothing over different 
lengths we derived a set of PDF's for different values of the 
linear rms fluctuation <jl- 

The results are displayed in Fig. []]. It is clear that the 
SC approximation provides excellent agreement with the nu- 
merical PDF even for ax 1, except for the large-5 tail 
of the distribution. These results are consistent with those 
of Protogeros, Melott, & Scherrer (1997), who used a less- 
accurate Lagrangian mapping which mimics the SC approx- 
imation and found good agreement with the PDF for power- 
law initial conditions. 



3 THE REDSHIFT DISTORTION OF THE PDF 
3.1 Preliminaries 

A particle at a position r in real space, with local velocity v, 
will be measured to lie at a position vu\ in redshift space, 
given by: 



•"(*) = r + {v r /H)r 



(22) 



where H is the Hubble parameter, and f points radially out- 
ward from the observer. For simplicity, we use the "infinite 
observer" approximation, so that the redshift direction is 
taken to lie parallel to one of the Cartesian axes, which we 
take to be the x axis. Then we have: 

r w = r + (v x /H)k. (23) 

It is convenient to make the substitution u = v/H, so that 



From mass conservation, we obtain the relation between r^ z ) , 
the value of r\ measured at a particular point in redshift 
space, and the corresponding rj at the same (Lagrangian) 
point in real space, where we can treat the redshift-space 
distortion as a Lagrangian mapping: 



n 



(25) 



V(z) l + (du x /dx) 

(see also Hui et al. 2000). Equation (^) is exact (for the 
infinite-observer case) and can be used, in principle, to cal- 
culate the redshift-distortion of the PDF. Such a calculation 
can be carried out for the skewness (Bouchet et al. 1995; 
Hivon et al. 1995; Scoccimarro et al. 1999), but there is no 
simple way to apply it to the full PDF. 

Hence, we make an approximation which is appropriate 
to the SC model. First we take du x /dx to be given by 

t - M 

where = V • u. This relation between du x /dx and is what 
would be expected in the limit where shear (and vorticity) 
are neglected, and our approximation has the advantage that 
the statistical properties of are well-known in tree-level 
perturbation theory. We get 

n 



n{z) 1 + (1/3)0 

which gives, for S, 

_ 8 - (1/3)0 
(z) 1 + (1/3)0 

We can use the continuity equation: 



0, 



(27) 



(28) 



(29) 



to express as a Lagrangian mapping of 8l- 



(30) 



e _ 1_ *? _ _ . 5 L dr) 

Hrj dt rj dSi, ' 

where fa = d In D/d In a, comes from applying the chain rule 
to the derivative of the linear growth factor: Sl = D(t)5o. 

3.2 Hierarchical moments 

Consider first the case of an unsmoothed density field. For 
small fluctuations all we need are the coefficients of the cor- 
responding Taylor expansion: 

\ " V k(z) xk 



^ k\ 

k 



We find for the first few orders 

= 1 + e 

v i{z.) = v-i + 2ev2 + 2e 2 



(31) 



(32) 



^3 + 3e^ 3 + 6e 2 (2^ 2 -l)+6e 3 



va(z) = v i + 4<^4 + 24e 2 (v 3 + v\ — 5^/2 + 1) 
+ 246 3 (3^ 2 -2) + e 4 , 

where e = /n/3. Note that in redshift space Vi( z ) 7^ 1, be- 
cause of the "Kaiser" effect. 

First consider the hierarchical amplitudes. For the red- 
shift variance we have: 
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Figure 2. The ratio of the redshift to real space variance or . /cr 2 
as a function of smoothing scale in SCDM simulations (points 
with error-bars) as compared to the "Kaiser" results (short- 
dashed lines) and the linear SC model (long-dashed lines). The 
dotted vertical line marks the scale where the linear variance is 
unity, i.e. ixj, = 1. Error bars are 1 — cr. 



(l + 6) 2 aI+0(4) 



(33) 



Note that this gives a reasonable approximation to the 
"Kaiser" effect in Eq. ([l]). For example, for fn — » 1 (e — > 
1/3), the expression above gives cr^/crf, = 16/9 ~ 1.78, 
compared with crf z) /a 2 L = 28/15 ~ 1.87 in Eq. These 
results are compared with (fully non-linear) SCDM simula- 
tions in Fig. El 



For the redshift-space skewness Ss(x) and kurtosis S, 



at leading order, are: 

3^2(z) S3 



S: 



3(z) 



T(z) 



(1 



2eS 3 + 6e 2 / 2 \ 
9 + 0(a L ) 



5 4(z), 



(34) 



5 4 ( Z ) 



S4. + 4eSj + 3e 2 (Si + 45 2 /9 + 325 3 /3 - 8) 
(1 + c) 4 

+ (l + e) 4 

where 6*3 = 3^2 and 6*4 = 12^| + 4!/3 are the corresponding 
real-space leading order results. In the limit fn — ► we 
get the real-space results, as expected. In the limit fn —* 1 
(e — > 1/3) we have, at leading order: 



S 3 « 

54( Z ) 



15 
16 



S 3 



(35) 



27 3 c2 27 9 
-S4 + -S3 + -S3--. 



So redshift distortions are not very large at leading order 
even in the case when fn — * 1. For the skewness, we get 
Smz) = 4.93, in excellent agreement with the results of exact 
perturbation theory, which predicts S3U) = 5.03 (Hivon et 
al. 1995), a difference of roughly 2%. For the kurtosis, we 
obtain 5*4( z ) = 47.5, (no exact calculation of £4(z) has been 
done thus far), which compares with a real-space value of 
S4 — 45.9. Like the skewness, the kurtosis is only slightly 
increased by redshift distortions. 



In practice, one needs to reach very large scales (R > 
100 /i _1 Mpc) to get to the regime where leading order is 
a good approximation. This is illustrated in Fig. ^, which 
shows how the ratio a^/a 2 estimated from SCDM Nbody 
simulations (presented in Section ^.2j ), only reaches the lin- 
ear regime at R > 100 ft _1 Mpc, even when the simulated 
box is L = 400 /i _1 Mpc and the scale of non-linearity is 
8 h^Mpc. 

We should therefore consider next to leading order cor- 
rections. For the variance we have: 



2 

O-(z) 



S 2,2(z) 
S 2,4(z) 



2 . 4 
s 2,2(z)T + *2,4(z)0"ii 



+ 

8 



(36) 



rS2,4 ' 



212S 2 + S 3 - 716 



For the skewness and fn 



162 
1: 



S3O) 

<S3,0(z) 
^3,2(z) 



S3,0(z) + Sz,2(z)VL, 

15 3 

16 S3 '° + 8' 

28 



(37) 



+ 



3772 - I6O8S3 



1149S| 



135*1 + 8645*4 



486 



where the values on the right-hand side are the ones in real 
space. We can see from the expressions above that the cor- 
rections are quite different in real and redshift space. 

Now consider the effects of smoothing. For power-law 
initial conditions, with a power spectrum P[k] oc k n , it is 
possible to show that (Juszkievicz, Bouchet and Colombi 
1993) 



53S = S3 + 7, 



(38) 



where 7 = — (n + 3). Bernardeau (1994) has shown that 
the above result is also true for a generic power spectrum, 
where 7 = f(R) = d\og(a 2 )/d\og(R). In Section 3.4 below, 
we will argue that the effects of smoothing on the redshift- 
distorted PDF are identical to the effects in real space (e.g., 
Eqs. pTj[ - ^l| ). Hence, we expect that 

S 3 (z)s = S 3 ( Z ) + 7 (39) 
Substituting these two equations into Eq. (^), we obtain 
53S + 2e53S+6e 2 + 7 6 2 ) 



S: 



3(z)S 



+ 



In the limit where fn — > 1, this becomes 



S 3 (z)s — Yq^ 3S ' 



3 - + ^L. 

8 16 



(40) 



(41) 



Since 7/ 16 is generally small compared to the other terms in 
the equation, smoothing has a minimal effect on the relation 
between the real and redshift space skewness. 

Fig. |§] shows a comparison of the leading order contri- 
bution to 53 in redshift (dashed line) and real (solid line) 
space as a function of smoothing scale R (which enters in 
the above expressions through 7 = 'y(R)). As can be seen in 
the figure, the real space predictions work well for R > 10, 
while in redshift space we need to go to larger scales to find a 
similar agreement with leading order results. (This is similar 
to what happens for the variance, e.g., see Fig. fy. 
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4 6 8 10 20 
R (h- 1 Mpc) 



4 



Figure 3. Comparison of the redshift (open squares) and real 
(closed square) space skewness S3, with 1 — a error bars, as a 
function of smoothing scale R, in SCDM simulations. Solid and 
dashed lines show the leading order SC result in real and redshift 
space. 



3.3 Dimensionless moments 

In discussions of large-scale structure, one is usually inter- 
ested in the hierarchical moments discussed in the previous 
section. However, as we shall see in Section 3.5 below, we find 
a simple empirical relationship between P[S] and P [<5(z)l- 
This result can be understood by considering the moments 
of the normalized field fi = 8/0. For the variance we have: 



(42) 



in both real and redshift space, by construction. We can 
define a dimensionless skewness: 



B3 = (// ) c = S3 cr, 
and in redshift space we find 
S3 + 2e5 3 + 6e 2 



B 



3(z) 



1 + e 
(1 + 2e) 



tl + O (crl) 



(43) 



(44) 



1 + e 

which in the limit /n 



B3 + 



_6£_ 

1 + e 
1 gives: 



o L + (at) , 



B 3M = \B 3 + ±a L + (at) 



(45) 



which should be compared with the leading order real-space 
result B3 = S 3 ol- At this order redshift distortions on the 
skewness are therefore still relatively small, while by defi- 
nition they are zero in the variance. For the dimensionless 
kurtosis we find: 

B Hz) = I B 4 + -L Bl + (35 3 - l)a 2 L + O (at) (46) 

The next to leading order term, B 3 ^a\ is, in the limit fn — > 
1: 



B 



3,3(z) 



2.745 3 + 4.25S3 + 4.73S'I 



2.95 

5.85 3 5 4 + 1.245 



I.6S4 



(47) 



while in real space (fn — > 0): 



B, 



-2 + 1 5 3 + | 51 

54 + To Ss 



ig 53 3 + y5 3 5 4 



(48) 



For the SCDM model at large scales, 
have: S3 ~ 2.86, 5 4 ~ 13.9 and 5 5 : 



B 



3,3(z) 



7.4 while in real space: B3 



where 7 ~ —2, we 
- 97. We then find 
~ 1.16. So there is 



a significant difference in these higher corrections, although 
the resulting change in B3 gets smaller as ox — > 0. 



3.4 The redshift-space PDF 

To obtain the final density at a given Lagrangian point, we 
first apply the SC approximation to obtain rj as a function 
of 8l in real space and then apply the redshift mapping in 
equation (p7|), with 8 given by equation (pcj), to go from 
rj to r](z). The combination of these two mappings can be 
treated as a single local Lagrangian mapping. We obtain, for 
the local Lagrangian relation between 5l and the redshift- 
distorted T)( z y. 



f(Sr 



Sl df(S L ) 



(49) 



where the function / is given by equations (^) and (|lo|). 
Finally, we must include the effects of smoothing. To do 
this, we apply equation ([TJ)]) or (|2l]) to the redshift-distorted 
density field. For a power-law power spectrum, for example 
we obtain: 



Sl(v M s/N) = F-'in^s/N^s/N) 



(n+3)/6 



while in the general case, we get 
Sl(v MS /N)=F- 1 (t 1(z)3 /N) 



<t(Ro) 



(50) 



(51) 



In these equations, F^ 1 is the inverse of the function F, 
given in equation ([49]), which transforms Sl into r\( z )- 

A couple of comments are in order regarding our 
smoothing methodology: we have convolved the local La- 
grangian mappings Sl — > V and rj — > rju) into a single local 
Lagrangian mapping r]^ = F{5l), given by equation (p7j[); 
it is this mapping which enters into Bernardeau's smoothing 
formalism. Note also that we are performing our operations 
in the correct order: equations ( [il^ ) and ( |50[ ) or (|5^) amount 
to first evolving the density field, then applying the redshift 
distortion, and smoothing the final result. 

For Gaussian initial conditions, our final expression for 
the redshift-distorted PDF is then: 



2-KOL 



■ exp 



2al 



d8 L 



V{z)S drj (z) c 



(52) 



where the function Sl = <5l(»7(z)s /N) is given by equations 
@ and @ or (|Tj), and N is given by 



N = 



2-KOL 



exp 



Sl(x) 2 



2al 



1 dS L (x) 
x dx 



dx. 



(53) 



We now test this approximation for the redshift- 
distorted PDF nu meri cally. The numerical simulations (de- 
scribed in Section 2/2) were used to calculate the redshift 
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space PDF for a variety of values of ox. These are com- 
pared with the SC prediction for the redshift-space PDF in 
Fig. g 

It is clear that the SC approximation gives a reason- 
able fit to the redshift-space PDF for ox ^ 0.4 but fails for 
larger values. Despite these problems, the SC approxima- 
tion appears applicable to larger values of ox than does the 
the Kaiser expression for (j( z ) (Eq. [j]). For example, Eq. (|l]) 
predicts a 2 z ^/a 2 = 1.87. In our numerical simulations, we 
obtain a 2 ^/a 2 = 1.21 for ox = 0.42, corresponding to the 
left panel in Fig. ^. Thus, Eq. ([[]) fails badly for gl as small 
as 0.4 (see also Fig. ^), for which the SC approximation 
provides reasonable agreement with the empirical PDF. 

3.5 An interesting empirical result 

For ox ~ 1, the SC approximation does a much better job of 
predicting the real-space PDF than the redshift-space PDF. 
Hence, we are motivated to consider other possible ways of 
predicting the latter. Note that when we normalize the dis- 
tributions to /i = S/a, the real and redshift space moments 
are very similar to leading order. The variances are identi- 
cal by construction and we have seen that the dimensionless 
skewness B$ and kurtosis B4 are quite similar at leading 
order. Thus one would expect that when at —> the dimen- 
sionless redshift PDF would be similar to that in redshift 
space. Note that fj, = 5/a and the dimensionless moments 



Table 1. The real-space and redshift-space values of the vari- 
ance, skewness, and kurtosis measured in SCDM simulations for 
a smoothing scale of R = 15 Mpc/h. 



a 


S3 


Si 


B 3 


Bi 


real 0.482 


2.66 


11.6 


1.28 


2.69 


z 0.510 


2.30 


6.8 


1.13 


1.77 



are the natural variables to quantify the amount of non- 
Gaussian effects on the PDF, as it is evident from the Edge- 
worth and Gamma expansions (e.g., see Gaztanaga, Fosalba 
& Elizalde 2000). 

In Fig. [| we provide a comparison of P[S / a]d[S / a] in 
real and redshift space. It is obvious from these figures that 
the agreement is excellent for ox < 0.6 and only breaks 
down at ox ~ 1- For ox < 0.6, the major difference occurs 
at the large-5 tails of the distributions, indicating that the 
higher order terms are becoming important. In particular, 
the real and redshift-space distributions diverge at S/a > 4. 

Note nevertheless that the agreement seems better than 
one would expect. Consider, for example, the values of the 
variance, skewness and kurtosis at R = 15 Mpc/h given in 
Table 1. We see that a^/a 2 lies below the linear theory 
prediction in Eq. ([[]) (or the SC prediction in Eq. |3oj ), while 
5*3(2) and S4( z ) are smaller than the corresponding real space 
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10 



-2 



4 

6/a 
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4 

(5/(7 
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Figure 5. A comparison of the real-space PDF, P[S/a] (solid curve) with the redshift-space PDF, P Is^/a^ j (dashed curve) 
SCDM simulations at the indicated values of a^. Error bars have been suppressed for clarity. 



values: in the opposite direction to the leading order predic- 
tion in Eq. ( [35| ) . These two effects seem to cancel each other 
so that the values of the dimensionless skewness Bz and 
kurtosis Bi in real and redshift space are quite similar. This 
suggests some additional cancellation of the redshift distor- 
tions which is not accounted for by the SC model, most likely 
due to nonlinear effects. 

Because of the potential significance of this result, we 
have repeated this calculation for a flat ACDM model with 
SIa = 0.8 (while this value for Qa is somewhat larger than 
that suggested by recent observations, this model and the 
SCDM model should serve to bracket the likely effect of red- 
shift distortions) . Our results are displayed in Fig. ^| Again, 
we see that, to a good approximation, P[S / a]d(S / a) is iden- 
tical in real and redshift space. The agreement here is even 
better than for the SCDM model; the two distributions are 
nearly identical for a as large as 0.8. 



4 CONCLUSIONS 

The spherical collapse (SC) approximation for the PDF 
can provide good agreement with both the real-space and 
redshift-space PDF from numerical SCDM simulations, but 
it fails at a much smaller value of <jl in the latter case. In real 
space, we get good agreement with the numerical results for 
ff L ~ 1) while in redshift space, good agreement can be ob- 
tained only for ox < 0.4. This is most likely due to nonlinear 



effects in the redshift-space PDF. An interesting advantage 
of the SC approach to derive the PDF is the possibility of 
extending the results to the case of non-Gaussian initial con- 
ditions (see also Gaztanaga & Fosalba 1998, Gaztanaga & 
Croft 1999). We only need to replace Po in equation (113) by 
the corresponding non-Gaussian initial PDF. 



We also find an unexpected result in our simulations: 
the redshift space PDF, P[<5( z j], is, to a good approxima- 
tion, a simple rescaling of the real space PDF, P[S], i.e., 
P[5/a]d[5/a] = P [<5( z )/<7 (z) ] d [<5 (z )/cr (z) ] . This rescaling 
provides an excellent approximation to the redshift-space 
PDF for o l < 1 in both the SCDM and ACDM models. 
We have not found an entirely satisfactory explanation for 
this. While it would be an obvious result from linear the- 
ory, its validity extends well beyond the linear regime. For 
example, this rescaling applies very well to the regime for 
which P[S] is highly non-Gaussian and the Kaiser prediction 
for (T( z ) from linear theory (Eq. [j]) breaks down. This result 
seems to stem from a cancellation between the linear effect 
given by Eq. (|2^) and the nonlinear redshift distortions. It 
is perhaps ironic that this simple-minded rescaling gives a 
much better prediction for the redshift-space PDF than our 
sophisticated perturbation-theory approximation. 
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Figure 6. A comparison of the real-space PDF, P[8/a] (solid curve) with the redshift-space PDF, P [5(je)/f(»)] (dashed curve) in flat 
ACDM simulations at the indicated values of a^. Error bars have been suppressed for clarity. 
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